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Abstract. We present the results of a study of charged pion production in 12 C + 12 C collisions at incident 
beam energies of 1A GeV and 2A GeV using the HADES spectrometer at GSI. The main emphasis of 
the HADES program is on the dielectron signal from the early phase of the collision. Here, however, we 
discuss the data with respect to the emission of charged hadrons, specifically the production of tt^ mesons, 
which are related to neutral pions representing a dominant contribution to the dielectron yield. We have 
performed the first large-angular range measurement of the distribution of 7T mesons for the 12 C + 12 C 
collision system covering a fairly large rapidity interval. The pion yields, transverse-mass and angular 
distributions are compared with calculations done within a transport model, as well as with existing data 
from other experiments. The anisotropy of pion production is systematically analyzed. 



PACS. 25.75.-q heavy-ion collisions - 25.75.Dw charged pion spectra 
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1 Introduction 

The investigation of nuclear matter at high temperature 
and high density is one of the major research topics in 
modern nuclear physics. Nucleus- nucleus collisions at rel- 
ativistic energies offer the unique possibility to create such 
highly excited nuclear matter in the laboratory [IIU.3,4, 
5,6. The study of particle production as function of beam 
energy, system size and centrality of the collisions has been 
instrumental in the past for understanding the approach 
of the system towards equilibrium and the generation of 
flow phenomena, as well as for gaining information about 
the nuclear equation of state. Collisions of the light 12 C 
+ 12 C system represent a link between the elementary 
nucleon-nucleon reactions and the heavy-ion collisions of 
large nuclei. Important physics issues in this context are 
the degree of thermalization achieved, the role of the mean 
field and collective motion in such a small system. 

In the few-GeV energy range, pions are the only abun- 
dantly produced mesons. In heavy-ion collisions their spec- 
tra and yields are affected by collective effects like ther- 
malization, directed and elliptic flow, as well as by possible 
modifications of the properties of the baryon resonances 
they stem from, in particular the A |7,8J. The subtle inter- 
play of the phenomena which change the characteristics of 
pion production with respect to nucleon-nucleon interac- 
tions is indeed a challenge to theoretical interpretations. 

Best suited for a description of all phases of the com- 
plex dynamics of heavy-ion reactions are transport mod- 
els, based on microscopic kinetic theory. The reaction is 
simulated as a network of multiple elementary collisions 
optionally embedded in a mean field equipped with momen- 
tum-dependent potentials. These assumptions arc then 
tested by comparing the experimental observables with 
the model predictions and allow to an understanding of 
the reaction dynamics. Transport models achieve a re- 
markable success in describing bulk properties of the in- 
teractions over a large energy and system size scale (cf. [TJ 
2,3,4,5,6]). At the same time, however, for special chan- 
nels, problems are met in reproducing precisely the ex- 
perimental data. For a recent comprehensive discussion of 
various differential pion observables and their comparison 
with model calculations in the region of 1A GeV see [9]. 

The High Acceptance DiElectron Spectrometer 
(HADES) [TU], in operation at the heavy- ion synchrotron 
SIS18 at GSI, Darmstadt, is designed for high-resolution 
and high-acceptance dielectron spectroscopy in hadron- 
hadron, hadron-nucleus, and nucleus-nucleus reactions at 
beam energies in the range from 1A GeV to 2A GeV. 
Being a charged-particle detector, it is of course also an 
efficient device for hadron detection. First results from 
HADES on dielectron production in 12 C + 12 C reactions 
have been presented in [TT1H2]. The interpretation of the 
dilepton spectra relies on a precise knowledge of the dif- 
ferential yields of neutral pions which are the source of the 
bulk of the detected dielectron pairs, namely via the ir° 
Dalitz and photon decays. In these analyses, the tt° yields 
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were inferred from the charged-pion yields measured by 
HADES in the same 12 C + 12 C data samples. 

In this paper we present detailed data on charged pi- 
ons obtained from 12 C + 12 C collisions at I A GeV and 
2A GeV. For the first time, large intervals of rapidity 
(p» ±0.8 in y/ybeam for 2A GeV) and of centre-of-mass 
angle (—0.7 < cos(9 cms ) < 0.7) are covered. Our results 
are compared to UrQMD transport-model predictions and 
experimental data from other experiments. 

2 Experiment 

HADES [TU] is a magnetic spectrometer designed as second- 
generation device for measurements of e + e~ pairs. The 
spectrometer, schematically depicted in Fig. [TJ is segmen- 
ted into six identical sectors that cover laboratory polar 
angles between 18 and 85 degrees. Its large azimuthal ac- 
ceptance covers between 65% and 90% of 2tt at small and 
large polar angles, respectively. The analysis of charged 
pions presented here is based on the same setup as used 
in the lepton analyses of rTTiri"2"] . A fast hadron-blind Ring 
Imaging CHerenkov counter (RICH) is used for electron 
and positron identification. In the pion analysis, the RICH 
information was not used. Four planes of Multi-wire Drift 
Chambers (MDC I - MDC IV), together with a supercon- 
ducting magnet, form the magnetic spectrometer for track 
reconstruction and momentum determination. In the re- 
gion behind the magnetic field, a set of electromagnetic 
Pre-Shower detectors (at polar angles 18°— 45°) [13] and a 
time-of-flight wall [2] are installed which form the META 
(Multiplicity and Electron Trigger Array). The time-of- 
flight detector wall is subdivided into 2 regions: TOF (at 
polar angles 45° — 85°), consisting of 384 scintillator slabs 
of varying length, which are read out at both ends, with 
a time-of-flight resolution of <r = 150 ps, and TOFINO 
(at polar angles 18° — 45°), consisting of 24 scintillator 
plates readout on one end, with a time-of-flight resolu- 
tion of cr = 450 ps. The TOFINO is placed directly in 
front of the Pre-Shower detector, which provides precise 
position measurement. The TOF/TOFINO detectors are 
also used for fast charged-particle multiplicity measure- 
ments. Together with the Pre-Shower detectors they pro- 
vide additional lepton/hadron discrimination power and 
track coordinate measurements with a spatial resolution 
in the range from 14 to 25 mm. 

A fast data acquisition system is used together with 
a two-level trigger scheme [T5"lll6j : (a) LVL1 is based on 
a fast determination of the charged-particle multiplicity 
(M c /j) in the time-of-flight detectors, (b) LVL2 is based 
on a real time identification of electron and positron can- 
didates. All LVL2-triggered events were written to tape, 
as well as a part of LVL1 (regardless of the LVL2 decision) 
events (typically 10%) for normalization purposes, hadron 
analysis and monitoring of the trigger performance. For 
the analysis presented here, only LVLl-triggered events 
were processed. 

In the very first HADES physics run, the detector was 
operated using only the following sub-systems: the RICH, 
the two inner MDC planes, and the META, i.e. the two 
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Fig. 1. Cut through two sectors of the HADES spectrometer. 
The magnet coils are projected onto the cut plane to visualize 
the toroidal magnetic field. The average distance between the 
target and the outermost detectors is about 210 cm. 

outer MDC planes were not operational and the position 
measurements from the META were used for tracking. In 
this mode, the collision system 12 C + 12 C at 2A GeV 
was studied with a beam intensity of heam — 10 6 parti- 
cles/sec impinging on a two- fold segmented carbon target 
with thickness of 2-2.5% interaction length. 1.67T0 7 LVL1 
triggered events with M c h > 4 were analyzed in this study. 
In the event reconstruction, the track segments measured 
in the two inner MDC planes were correlated with hits in 
the META. 

In the second data taking period the 12 C + 12 C system 
was studied at 1A GeV. Then, for the first time, a high- 
resolution tracking mode exploiting also the outer MDC 
planes was available. In this measurement, a carbon beam 
of 10 6 particles/sec was focused onto a carbon foil of 3.8% 
interaction length. 1.62 • 10 7 LVLl-triggered events with 
M c h > 4 were used in this analysis. 



3 Data Analysis 
3.1 Simulations 

Artificial 12 C + 12 C events were generated with the UrQMD 
(vl.3b) transport code [17lll8| . The detector response was 
simulated with the help of a GEANT 3.21 based pack- 
age [TTJ] including the geometry and characteristics of all 
HADES detectors. Then the same LVL1 trigger condition 



(Mch > 4) was applied. The resulting raw data were pro- 
cessed in exactly the same way as the real data and used 
for efficiency corrections as well as to estimate systematic 
errors. Details on the different procedures are given in the 
corresponding subsections. For 1A GeV we have analyzed 
2.14 • 10 7 (LVL1) UrQMD events, for 2A GeV 2.07 • 10 7 
events, i.e. numbers being comparable to the amount of 
analyzed real data. For this sample size, the statistical er- 
rors of the 7r yields in the region of interest are negligible. 
Systematic errors are not addressed by the UrQMD it- 
self, but may be elucidated by comparison with data. For 
the description of the apparently complex process of pion 
production in heavy-ion reactions in this energy range, the 
precision of transport models is not expected to be better 
than 20% g]. 

Additionally, we simulated the n meson emission using 
a simple Monte-Carlo event generator PLUTO [20], which 
assumes a thermal source modified by a polar angular dis- 
tribution. The simulation parameters used for PLUTO - 
inverse slopes and anisotropics - were derived from our 
measured data. The generated rapidity distribution has 
been used to extrapolate the pion yields outside our de- 
tector acceptance. Varying the input parameters of the 
generator within their errors provides estimates for the 
systematic error of the extrapolation. 

3.2 Momentum reconstruction 

When traversing the spectrometer, charged particles are 
deflected in the magnetic field and at the same time they 
leave "hits" in the MDCs and META detectors. From this 
information, together with the known magnetic field, their 
trajectories are reconstructed and their momenta are de- 
duced. 

Two different tracking methods have been developed 
and were both used in the present analysis (see [lOj for 
details). The first one is the "kickplane" algorithm which 
uses the position information delivered by the inner MDC 
chambers and the META system. In this case the momen- 
tum resolution a p /p, dominated by the limited position 
resolution of META, has been determined in simulations 
to be ~ 2.5% at a momentum of 150 MeV/c, with a lin- 
ear increase up to 16% at 1400 MeV/c, with a weak de- 
pendence on polar angle. The second method is based on 
a Runge-Kutta trajectory integration routine [3T] which 
uses the information from all four MDC planes, resulting 
in a resolution a p /p ~ 1.8% at 150 MeV/c and ~ 4.3% at 
1400 MeV/c. To obtain the results presented in this pa- 
per, the kickplane method has been applied to the 2 A GeV 
data, and both methods were used and compared for the 
1A GeV data. 



3.3 Particle identification 

Particle identification in the HADES data analysis (for de- 
tails see PUj) is based on Bayesian statistics [2"2"I|2"3"] . The 
method allows to evaluate the probability that the recon- 
structed track can be related to a certain particle species 
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Fig. 2. Velocity vs. charge-times-momentum of charged particles as seen by the HADES detector from 12 C + 12 C collisions at 
2A GeV (left). Projection onto the velocity axis of positively charged particles with momenta 350±5 MeV/c and 6 — 60°±5° 
(right). Fitted signal and background distributions are shown as lines. 



(e.g. proton, kaon, pion, electron, etc.). It combines sev- 
eral observables from various sub-detectors (e.g. time-of- 
flight, energy loss) via probability density functions (PDF) 
determined for each observable and for all possible par- 
ticle species. The probabilities for different mass assign- 
ments of any given track are calculated from the assumed 
abundances of the individual particle species and from the 
PDFs of all measured variables. The latter ones are ob- 
tained from simulations. If the assumed abundances differ 
significantly from the final results, the procedure is re- 
peated with updated input distributions. It converges typ- 
ically after one or two iterations. The performance of the 
method in terms of efficiency and purity is evaluated in 
detailed simulations and simultaneous comparisons with 
the real data. In our case, hadron identification has been 
performed using measured momenta and corresponding 
velocities computed by means of the time-of- flight. For 
more sophisticated analyses, like electron or rare-hadron 
identification, data from the RICH and PreShower detec- 
tors as well as the energy loss in META and MDCs can 
be used in addition. 

The method used for Particle IDcntification (PID) is 
illustrated in Fig.Ofor the case of particle velocity (right) 
deduced from the measured time-of-flight and track length 
( "velocity- vs- momentum" algorithm). Particles with dif- 
ferent mass occupy different regions in the velocity-vs- 
momentum distribution (left side); the pronounced ridges 
correspond to positively and negatively charged pions, 
protons and deuterons. The very rare kaons are not visi- 
ble in this representation. The Bayesian PID method re- 
quires the determination of the probability density func- 
tions for each particle species. In the case of the velocity- 
vs-momentum algorithm used here, the PDF is the prob- 
ability distribution of velocity. For each type of particle 
it has been determined in bins of momentum and polar 
angle. In those velocity distributions Gaussian fits were 
used to obtain the signal (i.e. particle yields) and a 2 nd - 
order polynomial fit to obtain the background (due to fake 
tracks). The fitted distributions were normalized to unity. 



Fig. [2] (right side) shows as an example such a fit for the 
momentum bin 350±5 MeV/c in the polar angle range 9 
= 60°±5°. 

Two quality parameters are used to characterize the 
performance of the method [24]: the PID efficiency and 
the PID purity. The PID efficiency et (p, 0) is the probabil- 
ity that a particle with the true type t is identified as type 
t. The PID purity n t (p, 9) is the probability that a parti- 
cle that is identified as type t is truly of type t. The PID 
efficiency and purity have been studied in detailed sim- 
ulations with events generated with the UrQMD model. 
The critical parameter here is the time resolution, which 
is however well known. Because of the moderate time res- 
olution of the presently installed TOFINO detectors this 
limits the region in which we can use the method for 7r + 
and p identification to momenta < 1000 MeV/c . Finally, 
we have checked that varying particle abundances by a 
factor as large as 2 does not change the results signifi- 
cantly. 



3.4 Total correction 

Fig. [3] shows the dependencies of the PID efficiency and 
purity on momentum for 7r + , 7r~ and protons for the 
2A GeV data in the TOFINO (right) and TOF (left) 
region. The efficiency of pion and proton identification 
is larger than 95% for all momenta in the TOF region. 
In the TOFINO region, with its reduced time resolution, 
the efficiency to identify positively charged pions drops 
steeply above 1000 MeV/c due to the ambiguity with the 
protons. The purity of pions (lower plots) does not reach 
unity mainly because part of the tracks identified as pions 
are muons from in-flight pion decays. The strong contam- 
ination of the positively charged pions with protons for 
momenta above 1000 MeV/c is again due to the low time 
resolution of the TOFINO. 

After the particle identification is done for all tracks, 
the resulting yields are corrected for efficiency and pu- 
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Fig. 3. Efficiency (top) and purity (bottom) of the PID method versus momentum for -k^ and protons in two sub-systems of the 
HADES detector: TOF (left) and TOFINO+PreShower (right) by using the kickplane reconstruction algorithm in 12 C + 12 C 
collisions at 2 A GeV. 



rity of the PID method, as well as for the detector and 
tracking efficiencies. The detection/tracking efficiency has 
also been obtained from Monte Carlo simulated and re- 
constructed UrQMD events. The reliability of the simula- 
tion was cross checked with real data by comparing and 
matching the information from different detectors and by 
selecting clean samples of (elastic) elementary collisions 
in which the track configurations can be obtained from 
the trigger detector (META) information only [10] . Using 
this information, realistic parameters were validated for 
the simulation. 

The total correction applied to the reconstructed par- 
ticle yields reads 



w t {p, 0) 



*t(p, 0) 



e t (p,6) ef et (p, 



(1) 



where ef et (p,6), the detection efficiency, subsumes detec- 
tor, track reconstruction and acceptance losses. It should 
be noted that we specify ef et (p,0) as function of and 
p, while averaging over the azimuthal angle. In this way, 
corrections for missing geometrical acceptance at some az- 
imuthal angles (namely the spaces occupied by the six 
magnet coils) are accounted for as well as the losses due 
to pion in-flight decays. 

Figure [4] shows the dependence of the detection effi- 
ciency for tt^~ and p as a function of momentum. The dif- 
ference between proton and tt efficiencies is again in part 
caused by the in-flight decay. 
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Fig. 4. The detection efficiency £t et (p,@) vs. momentum for 
■k^ 1 and protons using the kickplane reconstruction in 12 C+ 12 C 
collisions at 2A GeV. 



The total correction is applied to the data for each 
momentum and polar angle bin, and for each individual 
particle species. This is done only for bins with sufficiently 
high efficiency ef et {p,9) > 0.35 in order to avoid large 
corrections at the sector boundaries. Data outside of this 
fiducial volume were excluded from further analysis. 



6 G. Agakishiev et al.: Measurement of charged pions in 12 C + 12 C collisions at 1A GeV and 2A GeV with HADES 



C 

a> 0.8 
> 

0) 

"I 0.6 

u 

<X> 

§ 0.4 

U 
T3 

Z 0.2 
T3 



7l + 

D emitted 
• corrected 
° identified 



,ooooooooooo c 



u -1 -0.8-0.6-0.4-0.2 -0 0.2 0.4 0.6 0.8 1 

COS(9 cms ) 

Fig. 5. Ccnter-of-mass polar angle distributions of ix + . The 
squares represent the pions as generated by UrQMD. Only n 
mesons with p cms > 200 MeV/c have been selected. The open 
circles show those generated pions which are detected and iden- 
tified in the HADES acceptance. The full circles depict the 
result of the efficiency and purity correction to the accepted 
simulated pions. 

Figure [5] presents simulated polar distributions of pi- 
ons in the center-of-mass system (cms) . It shows the iden- 
tified 7r + before and after applying the total correction, 
together with the primordial distribution delivered by the 
UrQMD model for 2A GeV 12 C + 12 C collisions. This 
self-consistency check quantifies w t (p,0) as a function of 
cos 9 cms and demonstrates the wide coverage of our spec- 
trometer. In Fig.[5]the angular anisotropy of pion emission 
in the UrQMD model is clearly visible. 

The point-to-point systematic errors of the measured 
distributions stem mainly from imperfect modelling of the 
detectors, reconstruction and identification efficiency. They 
are estimated as « 5%, based on a comparison of measure- 
ments in the six independent HADES sectors. It should be 
noted that although the six sectors are identical in design, 
they are different at small scales like malfunction of elec- 
tronic channels and small mechanical displacements. The 
resulting errors are of statistical nature and thus appro- 
priate for fits to the various distributions. The geometry 
of the chambers and the magnetic field measurement have 
some influence on high resolution tracking, but do not 
affect continuum distributions as those discussed in this 
work. 



3.5 Event selection 

For the present analysis, we used the HADES LVLl-trig- 
gered events, which are characterized by a hit multiplicity 
M c h > 4 in the time-of-flight detectors. The correlation 
between the LVL1 trigger condition and the centrality of 
the reaction has been studied in Monte-Carlo simulations 
using the UrQMD and GEANT codes. Fig. [6] shows the 
simulated impact-parameter distributions. For the "mini- 
mum bias" events corresponding to the total reaction cross 



section we require in UrQMD at least one nuclear inter- 
action (distribution marked by circles in Fig. [S]). Then 
we pass these events through our analysis code and re- 
quire that they fulfill the LVL1 condition (triangles in 
Fig. [S]). We found that the LVLl-triggered events corre- 
spond to 52% and 60% of the total reaction cross sec- 
tion in 12 C + 12 C collisions at 1A GeV and 2 A GeV, re- 
spectively. To extract the average number of participants 
for our trigger-biased events we proceeded in the follow- 
ing way: For minimum-bias events the average number 
of participating nucleons was estimated from a geometri- 
cal model [25], namely (A™- rt ) = A/ 2 = 6. For reactions 
accepted by LVL1 we deduce the mean (A part ) by com- 
paring the pion multiplicity of UrQMD for LVLl-accepted 
events to the one of minimum-bias events, using {A part ) 
scaling of pion production (see e.g. figure 32 in [2j), i.e. 
(Ap art ) = 6- (M^ VL1 )/(M^ m - b -). The LVL1 trigger effect 
is significant, and the number of participants increases by 
~ 40%. The average impact parameters, the average pion 
multiplicities and average number of participants from 
UrQMD are listed in Table[T]for true minimum-bias events 
and after applying the LVL1 trigger. 



Table 1. Average impact parameters, pion multiplicities, and 
average number of participating nucleons from UrQMD calcu- 
lations for 12 C + 12 C at 1A GeV and 2 A GeV before and after 
applying the LVL1 trigger condition. 
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<6> 


{M«+) 


<M„-) 




minimum- 


bias events 3.95 fm 


0.83 


0.83 


6 


LVLl-trig 


gered 3.18 fm 


1.15 


1.17 


8.38 



The distributions of the number of reconstructed tracks 
per LVL1 event of data and UrQMD simulations are in 
reasonable agreement, as shown in Fig. [7] This is a con- 
firmation that the modelling of the detector and tracking 
efficiency as well as of the LVL1 event selection in our sim- 
ulation is realistic. From the differences of the measured 
and simulated distributions of the number of charged hits 
in META and the number of reconstructed tracks we es- 
timate a systematic error for the mean number of partici- 
pants determined above, resulting in (A part ) = 8.61 ±0.60 
and 8.38±o 17 for 1A GeV and 2A GeV, respectively. 

Comparing the pion momenta and angular distribu- 
tions for minimum-bias and LVL1 UrQMD events is an 
important check for how the trigger condition used in the 
experiment influences the measured spectra. We do not 
observe significant deviations between the two data sets, 
e.g. for the transverse mass distribution, the difference in 
the inverse slopes (see Section 4.1) is less than 1 MeV, 
and the Ai parameters of the polar angle distributions 
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(see Section 4.4) differ by less than 10%, while their mo- 
mentum dependencies arc identical. The reason is that for 
a small system such as 12 C + 12 C it is basically impossible 
to select events with well defined centrality by a charged- 
particlc multiplicity trigger. This fact is clearly seen in 
Fig. [6l which shows a very broad distribution in impact 
parameter from the LVL1 events. 



4 Results 



Figures [8] and [9] exhibit the measured and simulated 
transverse mass distributions of ir + and 7r~ in different 
intervals of normalized rapidity y = (yi ab - y cm s)/Vc ms 
for 12 C + 12 C at 1A GeV and 2 A GeV, respectively. The 
systematic errors of the data are estimated from point-to- 
point differences between distributions from the six inde- 
pendent HADES sectors as ~ 5% (see Sect. 3.5). 

The transverse- mass (mi) distributions have been fit- 
ted for each rapidity bin using one or two exponential 
functions. The fit with two slopes employs 



4.1 Transverse-mass distributions 1 dN(y) „ , . / m± \ „ , . / m± 

- C x {y) exp -— - +C 2 (y) exp - 



dm ±- V Ti{y)J V T 2{y)J 

(2) 
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Fig. 8. Transverse-mass distributions for positively (left) and negatively (right) charged pions in different slices of rapidity 
derived from the data in 12 C + 12 C collisions at 1A GeV (LVL1 "semicentral" events). Full lines show the results of fits to 
the data (symbols) using one exponential function, while dashed lines show fits of the UrQMD distributions using the same 
fit function. Error bars (systematic and statistical ones) are not visible at this scale. Both data and UrQMD distributions are 
normalized to the number of LVL1 events. 



with mj_ = (pl+m 2 ) 1 / 2 , and p± as transverse momen- 
tum. C\ t 2 are normalizations and Ti^ the inverse slope 
parameters; these parameters depend of course on rapid- 
ity. For the 2A GeV data, the two-component fit describes 
the experimental data better than a fit with one slope only 
(i.e., C2 = 0). For example, for the data at midrapidity 
X 2 /ndf is ~ 1.0 as compared to 4.8 for the fit with one 
slope. 

Figure clearly demonstrates that for the lower bom- 
barding energy of 1A GeV, a fit with one slope is suffi- 
cient for the description of the spectral shape. The inverse- 
slope parameters for n mesons at midrapidity for 1A GeV 
(-0.15 <y < 0.15) and 2A GeV (-0.1 < y < 0.1) are 
summarized in Table [2] using one or two exponential func- 
tions. The slopes of 7r + and tt~ agree within error bars for 
the single exponential fit. At 2A GeV, UrQMD predicts 
different spectral shapes (purely exponential as compared 
to the concave shape of the real data), while at 1A GeV 
the agreement of UrQMD with data is better. Note that 
the errors of the fitted parameters are much larger in the 
case of two-slope fits than for the single exponential fits. 
The reason is occurrence of large correlations between the 



fit parameters in the two-slope case, which cause a con- 
siderable increase of the errors. 

Thanks to the wide acceptance of the present experi- 
ment we can directly compare the data with results of pre- 
vious experiments on charged pion production in 12 C + 12 C 
collisions at 1A GeV and 2A GeV done by the KaoS collab- 
oration [26] , and for neutral pions at 1 A GeV and 2 A GeV 
obtained by the TAPS collaboration [57] ■ For this purpose 
our data were passed through the respective acceptance fil- 
ter of the previous measurements: ±4° around the labora- 
tory polar angles 9i a b of the KaoS setup, and 0.42 < yi a b < 
0.74 and 0.80 < y lab < 1.08 for 1A GeV and 2A GeV, re- 
spectively, for the TAPS case. Then we rescaled the mul- 
tiplicity per one LVL1 event to the minimum-bias cross 
section. For this we used a total reaction cross section of 
0.95 b (calculated according to irr 2 (A^ 3 + A 1 J 3 ) 2 , assum- 
ing ?*o = 1.20 fm) and the ratio of (A part ) for minimum- 
bias and LVL1 events from Table Q] Fig. ITD1 shows a com- 
parison of our pion transverse-mass and momentum dis- 
tributions with those of the TAPS and KaoS experiments, 
respectively. The TAPS data are shown both in linear and 
logarithmic scale to facilitate the comparison of yields as 
well as of spectral shapes. 
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Fig. 9. Transverse-mass distributions for positively (left) and negatively (right) charged pions in different slices of rapidity 
derived from the data in 12 C + 12 C collisions at 2A GeV (LVL1 "semicentral" events). Full lines show the results of fits to the 
data (symbols) using two exponential functions, while dashed lines show fits of the UrQMD distributions using one exponential 
function. Error bars (systematic and statistical ones) are not visible at this scale. Both data and UrQMD distributions are 
normalized to the number of LVL1 events. 



It is apparent that the yields and spectral shapes mea- 
sured by all three experiments arc in general fairly similar. 
The differences between integrated yields are within errors 
and do not exceed 10%. Looking at Fig.Qjllin more detail, 
one sees that the measured charged-pion data agree for 
both bombarding energies (lower two plots, see insets for 
data ratios) over a large momentum range within errors 
(dominated by large systematic errors of the KaoS data) . 
On the contrary, the comparison of charged and neutral 
7T mesons (upper 4 plots) shows significant deviations for 
both beam energies, most apparent in the 2A GeV case 
(see inset of the upper right plot). As a consequence, our 
data on charged pions can be described only by a two- 
slope fit, while the TAPS ir° data at this energy need only 
a single-exponential fit [27] • This may indicate differences 
in the reaction dynamics of charged and neutral 7r mesons, 
not present in transport codes. It is interesting to note that 
the pattern of the 7r°/7r~ ratio changes with beam energy 
as well (see insets in middle row in Fig. [TO]) . 

Other fits of our pion transverse-momentum distribu- 
tions arc in principle conceivable. For instance, a blast 
wave fit (see e.g. [28]) can be used. However, if applied to 



only one particle type it does not allow for an unambigu- 
ous determination of the flow parameter. 



4.2 Rapidity distributions 

As seen in the previous section, the HADES acceptance 
in transverse momentum is rather large. For the missing 
parts of the acceptance at low and high p± we extrapo- 
lated the pion yield. In doing so, for each slice centered 
at rapidity y the corresponding p± distribution was fitted 
with the function 

= Cl(y) exp i~m) + C2iy) cxp (rm) ■ 

(3) 

The fit results were then used to integrate the tt ± yield 
outside the acceptance. The resulting correction was » 5% 
and « 2% for ir + and ir~ respectively, except for data at 
large rapidities, where the correction was 10 - 20%. Fig.fTTI 
shows the rapidity distributions normalized to the num- 
ber of LVL1 events obtained from the integration of the 
extrapolated p± spectra of tt + and tt~ from 12 C + 12 C 
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Fig. 10. Comparison of transverse-mass and transverse-momentum distributions of it mesons measured in the present experi- 
ment (HADES) and the previous TAPS [27] and KaoS |26] experiments. Top: Transverse-momentum (p±) distributions of ir ± '° 
mesons for 1A GeV (left) and 2A GeV (right). Center: Transverse-mass (mi) distributions of n ± '° mesons for 1A GeV (left) 
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GeV (left) and 2A GeV (right). The HADES data were rescaled to minimum- bias cross sections and filtered with the KaoS and 
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at 1A GeV and 2A GeV. UrQMD simulation data (again 
normalized to the number of LVL1 events) are displayed as 
well. In order to extrapolate our data outside the accepted 
rapidity range, we used the Monte-Carlo event generator 
PLUTO described above. The simulation parameters - in- 
verse slopes and anisotropies - were derived from our data, 
and the resulting rapidity distributions were normalized to 
the measured data. 



The rapidity distributions exhibit a Gaussian-like shape 
with a standard deviation (<r referring to the scaled rapid- 
ity yo) of about 1.0 for both beam energies. In the case of 
the 1A GeV data, the experimental rapidity distribution 
is about 20% narrower than the one from UrQMD. This 
is in agreement with the finding (see section 4.4 below) 
that the anisotropy parameter has a lower value in the 
experimental data than in UrQMD simulation. Note that 
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Table 2. Inverse slope parameters in units of MeV for it measured at midrapidity derived from the data (using one and two 
exponential functions) and UrQMD (using one exponential function) in 12 C + 12 C collisions at 1 and 2A GeV. 







Beam 


energy = 1A GeV 






Particle 


Data 
T (1 slope) 


X 2 /ndf 


UrQMD 
T (1 slope) X 2 /ndf 






7V+ 


57.8 ± 0.3 


1.7 


55.4 ± 0.3 2.2 








57.9 ± 0.3 


1.4 


55.4 ± 0.3 2.0 










Beam 


energy = 2A GeV 






Particle 




Data 


UrQMD 




T (2 slopes) 


X 2 /ndf 


T (1 slope) X 2 /ndf 


T (1 slope) 


X 2 /ndf 


7V+ 


47.7 ± 6.2; 
90.6 ± 3.3 


0.9 


80.9 ± 0.5 4.7 


86.5 ± 0.6 


1.5 


7T~ 


46.4 ± 5.2; 
84.4 ± 2.1 


1.2 


76.7 ± 0.5 4.9 


86.7 ± 0.6 


1.4 



for the 2A GeV slight underestimation of our data 

by UrQMD is observed. 



4.3 Multiplicities 



Table 3. Particle yields per reaction (LVL1 trigger condition) 
of vr* from 12 C + 12 C collisions. N„ and N^(4tt) are the mea- 
sured and 4ir extrapolated yields, respectively. The statistical 
errors are negligible. Shown are the systematic errors due to 
the various efficiency corrections and to the 4n extrapolation 
(see text for details). 



beam energy 


particle 


N n 


N^{4n) 


(A GeV) 








1 


7T+ 


0.36±0.02 


0.46±0.03±0.05 


1 


TT~ 


0.38±0.02 


0.49±0.03±0.05 


2 


7T+ 


0.77±0.04 


1.19±0.06±0.11 


2 


TT~ 


0.82±0.04 


1.28±0.06±0.12 



Pion yields (N^) per reaction (under the LVL1 trigger 
condition) within the HADES acceptance region as shown 
in Fig. [TT1 and in full phase space are presented in Tabled 
The systematic error of the measured yield due to uncer- 
tainties in the detection/reconstruction/identification ef- 
ficiency is estimated as 5%, based again on a comparison 
of measurements in the six independent HADES sectors. 
The extrapolation of the yields to full phase space is based 
on the integration of the rapidity distributions simulated 
with PLUTO (see above) and normalized to the data in 
the rapidity range covered by HADES (see histograms in 
Fig. [TTj) . Varying the input parameters of the PLUTO 
simulations - inverse slopes and angular anisotropy pa- 
rameters - within their experimental errors, the differences 
between the rapidity distributions give us an estimate of 
the systematic error of the yield extrapolations of 9%. 



Using the estimated averaged number of participants 
in the LVLl-triggered events (see in Table [1]), the result- 
ing 7r~^~ multiplicity per participant (averaged for tt + and 
7I-) is then 0.055 ± 0.003 ± 0.005 ± 0.004 and 0.147 ± 
0.007 ± 0.013 t° , Q21 at 1 A GeV and 2A GeV, respectively. 
The three systematic errors of the yields correspond to 
uncertainties connected with the efficiency/purity correc- 
tions (5%), the extrapolation to full solid angle and full 
kinematic phase space (9%), and the determination of the 
number of participating nucleons (see Section 3.5). 

Table [4] shows a comparison of measured 7r meson mul- 
tiplicities per participant and results of our UrQMD sim- 
ulations as well. The results agree within quoted errors. 

Table 4. Comparison of multiplicities per participant of tv 
mesons derived from our data with UrQMD results. 



beam energy 


particle 


this work 


UrQMD 


(A GeV) 




xlO -3 


xlO -3 


1 


l/2(7r++7r-) 


55 ±3 ±5 ±4 


59 


1 


TT 




67 


2 




147 ± 7 ± 13 ±21 


137 


2 


7V° 




159 



As we have shown above by comparing the momentum 
distributions, our results on integrated charged pion yields 
are consistent with the TAPS 7T° as well as the KaoS ir + 
data. Due to experimental errors, it is however difficult to 
draw a conclusion on the difference in production yields of 
neutral and charged 7r mesons predicted by UrQMD (see 
Table SJ). 

4.4 Angular distributions 

The measured centre-of-mass polar angular distributions 
of pions produced in 1A GeV and 2 A GeV 12 C -I- 12 C col- 
lisions are exhibited in Fig. together with the corre- 
sponding UrQMD distribution. Pions with centre-of-mass 
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Fig. 11. The rapidity distribution of positively (left) and negatively (right) charged pions produced in 12 C + 12 C collisions at 
1A GeV (top) and 2A GeV (bottom) for the LVL1 (semicentral) events. Circles with error bars show data, while full squares 
depict UrQMD calculations. The distributions obtained with the PLUTO generator are shown as histograms normalized to the 
measured yields (with ir + and ix~ averaged). In the 1A GeV case, the data points reflected from forward rapidities are shown 
as open symbols. The slight asymmetries of dN/dy with respect to reflection at yo = are used to check the systematic errors. 



momenta between 200 and 800 MeV/c have been selected. 
No losses in acceptance occur in this phase space region 
for the range of polar angles shown in Fig. [T2j The sys- 
tematic errors of the data are again estimated from the 
point-to-point differences between distributions from the 
six independent HADES sectors to be w 5% (see Sect. 
3.5). As seen from the comparison of the measured data 
points and the ones reflected around 90° (see Fig. ITS")) , this 
error underestimates the observed differences between the 
forward and backward hemisphere by a factor of < 1.5. 



In the symmetric collision system 12 C + 12 C the polar 
distributions in the center-of-mass system can be fitted 
with the following expression 



dN 



d(cOS 9 cms) 



A 1 (l+A 2 cos 2 9 cms ). 



(4) 



The fit parameter A 2 characterizes the anisotropy of the 
pion source, and A\ is a normalization. As visible in Fig. 1121 



the data show strong anisotropies quantified by A 2 = 
0.88 ± 0.12 and 1.19 ± 0.16 for beam energies of 1A GeV 
and 2A GeV, respectively. The UrQMD model gives sim- 
ilar anisotropies with A 2 — 1.45 and A 2 = 1.12 (A 2 = 
0.56 and A 2 = 0.70 when integrated over all momenta, 
including also the region outside our acceptance). 

We observe a strong dependence of the anisotropy on 
momentum. This is evident from Fig. 1 1 31 where A 2 is dis- 
played as a function of the pion's centre-of-mass momen- 
tum from fits of Eq. ([4} to our data. The data points 
for 7t + at highest displayed momenta are not shown, as 
they correspond to the phase space region (p; & approach- 
ing 1000 MeV/c) where reliable particle identification be- 
comes gradually impossible (see Fig. |3]in Sect. 3.3). 



One can see that the anisotropy steadily increases with 
momentum for both pion charges and both beam energies 
up to A 2 ~ 1.0 — 1.5, around 400 MeV/c, where it has 
a tendency to level off. This behavior is fairly well re- 
produced by the UrQMD model, which, however, tends to 
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Fig. 12. Polar angle distribution in the center-of-mass system of positively (left) and negatively (right) charged n mesons 
produced in 12 C + 12 C collisions at 1A GeV (top) and 2A GeV (bottom) for the LVL1 (semicentral) events. Pions with center- 
of-mass momenta 200 - 800 MeV/c have been selected. Full circles show measured data. The data points reflected around 90° 
are shown as open squares. The full lines show the fit by Eq. (|4} , the extrapolations of the fits outside the acceptance are plotted 
as dashed lines. The UrQMD points are displayed as triangles, no fits to the UrQMD points are shown. 



level off at somewhat larger values of the anisotropy in case 
of the 1A GeV data (see Fig. [T3|) . For the data in the re- 
gion 100 MeV/c < p cm s < 200 MeV/c, the anisotropy can 
be fitted in a limited O cms range only. The results show 
significantly lower A2 parameters for both the data and 
UrQMD. The close-to-zero (at 2 A GeV) and even slightly 
negative anisotropies given by UrQMD for low momenta 
at 1A GeV are not seen in the experiment. 

It is interesting to compare these distributions to the 
corresponding ones of the elementary NN — + N A — * 
NNir reactions which are expected to be the dominant 
source of pion production at these beam energies [HIE]. 
In order to estimate the latter we have again used our 
PLUTO generator employing measured A distributions to 
model the proton-proton reaction [29] . We find that the 
shape of the simulated A2 distributions is similar to the 
one found in 12 C + 12 C, but A2 levels off at substantially 



higher values: namely 3.5 and 5 for 1A GeV and 2 A GeV, 
respectively. The lower A2 values found for 12 C + 12 C 
may be attributed to pion re-scattering and final state in- 
teractions even in the small system under consideration 
here. 

Prior to the present work no data had been published 
on pion anisotropies in 12 C + 12 C collisions. Early studies 
of pion production in Ne-induced reactions at 0.8A GeV 
on NaF, Cu and Pb targets had found almost isotropic 
angular distributions for very low pion cms kinetic en- 
ergies (E n + < 50 MeV), but substantial anisotropies for 
all higher energies (E % + > 150 MeV) [30]. The two clos- 
est systems studied most comprehensively in this respect 
are Ar+KCl at 0.8A GeV as well as 1.8A GeV [3H[32] 
and Ca+Ca at 1.93 A GeV [9]. While in both cases simi- 
lar momentum-averaged anisotropies were observed, with 
values of (A2) = 0.5 — 0.6, only the Ar+KCl data display 
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Fig. 13. Dependence of the anisotropy parameter Ai on the momentum in the center-of-mass system for positively (left) 
and negatively (right) charged pions produced in 12 C + 12 C collisions at 1A GeV (top) and 2A GeV (bottom) for the LVL1 
(semicentral) events. Circles with error bars show fits to the data, while squares exhibit fits to UrQMD simulations after they 
were subjected to the full analysis chain. Errors of the UrQMD points are smaller than the symbol size. 



a strong pion-energy dependence of A2, peaking around 
E n = 200 - 300 MeV. Based on a comparison with trans- 
port calculations, the authors of Ref. [5] attributed these 
differences to the very different centralities covered by the 
two measurements. As seen from Fig. [13] and as discussed 
above, we do not observe in 12 C + 12 C, at both bom- 
barding energies, a rise and fall of A2, but rather a simple 
lcvclling-off with increasing pion cms momenta. 

As the momentum dependence of A2 has been claimed 
to depend strongly on the reaction centrality [9] , we tried 
to study this effect in our data sample. First we checked, 
how the distribution differs in UrQMD events for the mini- 
mum-bias and LVLl-triggered events (see Sect. 3.4). The 
effect is negligible, as could be expected for a very weak 
preference of more central events of the LVL1 trigger in 
the light 12 C + 12 C system. Then we selected the 10% of 
"most central" events in the measured and UrQMD data 
samples by requiring more than 3 identified protons in the 



event. Again, we did not observe any significant difference 
with respect to the LVL1 data samples. 



5 Summary 

In summary, the charged-pion characteristics in the reac- 
tion 12 C + 12 C at 1A GeV and 2 A GeV have been mea- 
sured in detail with the HADES spectrometer. The found 
7r meson yields are in good agreement with previous results 
obtained with the TAPS and KaoS detectors. The much 
larger acceptance of the present experiment, as compared 
to older measurements allowed for a more precise and re- 
liable extrapolation of the yields to full solid angle, which 
is essential for the normalization of our dielectron yields 
in this reaction system. 

Our data on the pion transverse-mass distributions at 
midrapidity can be described by a Maxwell-Boltzmann 
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function for 1A GeV, while the 2A GeV data show a strong 
second exponential component with smaller slope. This 
finding is in agreement with former results for the same 
system and energy for ir^ 1 [26], whereas n° data [27] ex- 
hibit only a one-slope behaviour. The reasonable agree- 
ment of our transverse-momentum spectra at 1A GeV 
with UrQMD calculations indicates that the degree of 
thermalization in the light 12 C + 12 C system is adequately 
reproduced in this transport model. At 2A GeV, however, 
some deviations of our data from the UrQMD results be- 
come visible. 

In contrast to the TAPS and KaoS data, which were 
measured in a limited angular range, strong pion anisotro- 
pics could be observed in the much larger acceptance re- 
gion of the present experiment. The systematics [9] of pion 
production in heavier systems at comparable beam ener- 
gies points to similar anisotropy values as extracted from 
our data sample. The asymmetries have a non-trivial mo- 
mentum dependence. A peaking of the asymmetry as re- 
ported in [31], however, can not be confirmed. The com- 
parison of our results with proton-proton data suggests 
that the angular anisotropies observed in 12 C + 12 C are 
remnants of the much stronger effect characteristic for in- 
elastic nucleon-nucleon scattering. 
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